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Using Monte Carlo techniques, the surface critical behaviour of three-dimensional semi-infinite 
ANNNI models with different surface orientations with respect to the axis of competing interactions 
is investigated. Special attention is thereby paid to the surface criticality at the bulk uniaxial 
Lifshitz point encountered in this model. The presented Monte Carlo results show that the mean- 
field description of semi-infinite ANNNI models is qualitatively correct. Lifshitz point surface critical 
exponents at the ordinary transition are found to depend on the surface orientation. At the special 
transition point, however, no clear dependency of the critical exponents on the surface orientation is 
revealed. The values of the surface critical exponents presented in this study are the first estimates 
available beyond mean-field theory. 



I. INTRODUCTION 

Critical phenomena at surfaces have been extensivjelK 
studied theoretically during the last three decades .EtEI 
The surface phase diagram of the three-dimensional semi- 
infinite Ising model with only nearest neighbor couplings 
is well established. Introducing two different ferromag- 
netic interactions depending on whether the neighboring 
spins are both located at the surface, Js > 0, or not, 
Jt, > 0, two typical scenarios are encountered. If the ra- 
tio of the surface coupling Jg to the bulk coupling Jb, 
r = Js/Jb, is smaller than pSu critical value, r^p « 1.50 
for the simple cubic latticeJj'Q the system undergoes at 
the bulk critical temperature Tc an ordinary transition, 
with the bulk and surface ordering occurring at the same 
temperature. Beyond this critical value one first observes 
the surface transition where the surface alone orders at 
a temperature Ts > Tc, followed by the extraordinary 
transition of the bulk at Tc- At the critical ratio r^p, the 
special transition point is encountered, displaying criti- 
cal properties distinct from those of the ordinary or the 
surface transition. The critical exponents of the differ- 
ent surface uriiversality classes are known with rather 
high precision .eI However, far less is known on the surface 
critical behaviour of models with competing interactions 
which is the subject of the present paper. 

Thfi-|p,axial next-nearest-neighbor Ising (ANNNI) 
modcH'tl is the best known of these models, its bulk prop- 
erties being intensively investigated since many years. 
Here, competition between ferromagnetic nearest neigh- 
bor and antiferromagnetic next-nearest neighbor cou- 
plings takes place in one direction. The Hamiltonian of 
the three-dimensional version, defined on a cubic lattice, 
may then be written in the form 



— ^Jb Sxyz {s(x + l)yz + Sx(y+l)z + Sxy(z + 1)) 



xyz 

-Kf) Jb ^ ^ ^xyz'^ 
xyz 



zSxy{z+2) 



(1) 



direction of competing interactions (here, the z-direction) 
is also called axial direction. 

Spatially modulated phases due to competing interac- 
tions are observed, anipiig others, in magnetic systems, 
alloys or ferroelectricsJ3^E3 Some compounds (as, for ex- 
ample, BCCDEj or NaV205EJ) display very rich phase 
diagrams with a multitude of commensurately or incom- 
mensurately modulated phases. A further remarkable 
feature of these systems is the possible existence of a 
special multicritical point called Lifshitz point .113 At a 
Lifshitz point, a disordered, a uniformly ordered and 
a periodically ordered phase become indistinguishable. 
Various, syjstems have been shown to possess a Lifshitz 
point Jia^lla^y 

Recently, there has been renewed interest in the uni- 
axial Lifshitz point encountered in thfi-piiase diagram of 
the three-dimensional ANNNI model.E3cZl This strongly 
anisotropic equilibrium critical point, located at si = 
0.270 ± 0.004 and ksT^/Jb = 3.7475 ± 0.0005Hj is 
characterized by the anisotropy exponent 6 = v^/v^ 
where and are the exponents of the bulk cor- 
relation lengths parallel and perpendicular to the z- 
axis. The value 9 « lO-iM has been obtained in a 
second-order e-expansionSo Recent thorough determi- 
nations of Lifshitz point critical exponents-uaing different 
techniques (fielitheoretical calculationtcjo and Monte 
Carlo methodscZl) have yielded excellent agreement. In 
addition, the scaling of the spin-spin correlator at the 
Lifshitz point was determinedcZI and its form shown to 
agree with the prediction resulting from a generalization 
of conformal invariance to. the strongly anisotropic scal- 
ing at the Lifshitz point .£3 

Whereas the properties of the bulk ANNNI model are 
well established, far less is known about the behaviour 



of the ANNNI mode: 
Monte Carlo studiesE 



he presence of surfaces. Early 
investigated ANNNI samples 



where Jt > and Kf, > are coupling constants. The 



with free surfaces but did not pay special attention to sur- 
face properties as they were exclusively interested in bulk 
behaviour. Thin ANNNI films with free surfaces were 
studied recently and they were shown to present a dis- 
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tinct phase diagram for every film thickness£2l An early 
attempt to study surface critical behawiour near the Lif- 
shitz point was undertaken by Gumba23 in the framework 
of mean-field theory. A much more elaborated mean-field 
treatment of the semi-infinite ANNNI modeL-sjps pre- 
sented recently by Binder, Frisch, and Kimball. L^tfl They 
considered two different surface orientations with respect 
to the axis ofj-cpmpeting interactions: surfaces oriented 
perpendicular^^ or parallels to this axis. For both cases 
the mean-field surface critical exponents at the Lifshitz 
points were determined and two different sets of critical 
exponents were obtained. This dependency of the critical 
behaviour on the surface orientation may be explained by 
the anisotropic scaling at the Lifshitz point. These au- 
thors also predicted the existence of a surface transition 
at the Lifshitz point-where only the surface orders, in 
variance with GumbsLJ who claimed that at the Lifshitz 
point the surface could not order before the bulk, ft 
is worth noting that no results beyond mean-field the- 
ory are available on the surface critical behaviour of the 
ANNNI model close to the Lifshitz point. 

In this paper I present the first Monte Carlo study 
of the surface critical behaviour at the uniaxial Lifshitz 
point encountered in the ANNNI model. This study uses 
a cluster flip algorithm, especially designed for simulating 
models with competing interactions, which has been suc- 
cessfully employed in a recent investigation of the bulk 
ANNNI modeLEZl 

The paper is organized in the following way. The next 
Section is devoted to the presentation of semi-infinite 
ANNNI models with two different surface orientations. 
Some details of the numerical method are also presented. 
Section III deals with surfaces perpendicular to the axis 
of competing interactions, paying special attention to the 
surface phase diagram and to the determination of sur- 
face critical exponents. Surfaces oriented parallel to this 
special direction are treated in Section IV. A short sum- 
mary and outlook conclude the paper. 

II. MODELS AND METHOD 

Due to the anisotropy of the ANNNI model, surfaces 
with different orientations are not equivalent. The two 
surface orientations with respect to the axial direction 
considered in this paper are the foUowings (see Figure 1): 
surfaces perpendicular to the axis of competing interac- 
tions (case A) and surfaces parallel to this axis (case B). 
The saipei-Oirientations were treated in recent mean-field 
studies.m 

For case A modified surface couplings with strength 
Js connecting neighboring surface spins are introduced 
in addition to the usual ANNNI interactions, see Eq. (1) 
and Figure la. Three different scenarios have to be dis- 
tinguished, depending on the value of the bulk compe- 
tition parameter Kf,. When Kf, is smaller than the Lif- 
shitz point value k^, the bulk undergoes a second order 



phase transition between the disordered high tempera- 
ture phase and the ordered, ferromagnetic, low tempera- 
ture phase at the critical temperature Tc{Kb)- This phase 
transition belongs to the universality class of the 3D Ising 
model. Consequently, the surface phase diagram will re- 
semble that of the 3D semi-infinite Ising model, with 
a possible shift of the location p£ the special transition 
point, rsp{Ki,), as function of k^O At the Lifshitz point, 
Kb = , the recent mean- field treatment yields a sur- 
face phase diagram for the semi-infinite ANNNI model 
similar to the Ising model, but with critical exponents 
which differ fromj-those of the Ising model. The con- 
jectured existenceE3 of a surface transition at the Lif- 
shitz point, and therefore of a special transition point, is 
in agreement with the Monte Carlo results presented in 
the next Section. Finally, for strong axial next-nearest 
neighbor bulk couplings, Kb > , the bulk phase tran- 
sition belongs to the universality class of the 3D XY 
model. Therefore, for weak surface couplings, the ordi- 
nary transition critical behaviour should be identical to. 
that of the three-dimensional semi-infinite XY model,L£l 
whereas for strong surface couplings the surface effec- 
tively decouples from the bulk and a two-dimensional 
bulk Ising critical behaviour is expected at the surface 
transition. These two critical lines meet at a multicritical 
point which should have interesting properties. However, 
I will in this study not consider the latter case as I am 
mainly interested in the surface critical behaviour at the 
Lifshitz point. 

Case B with surfaces oriented parallel to the axial di- 
rection (see Figure lb) is the most interesting but also 
the most demanding case. Introducing modified nearest 
neighbor, J^, and axial next-nearest neighbor couplings, 
Ks > 0, ya the surface layer leads to intriguing and very 
complex situations. For example, a multicritical point 
shows up where an ordinary transition with a modulated 
bulk meets a surface transition to a floating incommen- 
surate phase in the surface layer. In this study I will not 
discuss these pecularities but concentrate exclusively on 
the case Kg < nj^ . For these values of the surface compe- 
tition parameter Ks the corresponding two-dimensional 
ANNNI model presents a transition, belonging to the 2D 
Ising universality class, from the disorjdered phase to the 
low temperature ferromagnetic phase.E3 

Surface quantities are determined by simulating 
ANNNI films consisting of x Ly x spins, with free 
boundary conditions perpendicular to the surfaces and 
periodic boundary conditions otherwise. Hence, for case 
A the system consists of Lz layers with x Ly spins 
per layer, whereas for case B the samples are formed by 
Ly layers containing x Lz spins. The semi-infinite 
models are then obtained in the limit Lx, Ly, Lz — *■ oo. 
Some care is needed when choosing the shape of the fi- 
nite ANNNI films in the vicinity of the Lifshitz point. On 
the one hand, the special finite-size_effects coming from 
the ani^ropic scaling at this pointEZi are best taken into 
accounttZI by choosing anisotropic samples with an in- 
creased number of sites perpendicular to the axial direc- 



2 



tion. On the other hand, computation of surface quanti- 
ties is usually done by simulating samples where the lie-, 
ear dimension of the surfaces exceeds the film thickness .Ea 
This does not pose any problem for case A, as here the 
axial direction coincides with the direction perpendicu- 
lar to the surfaces. Therefore, for this case, systems of 
anisotropic shape with Lz layers and spins per layer 
are simulated, ranging from 10 to 120 and from 
10 to 60. For case B, however, surfaces are parallel to 
the axial direction. In order to balance the competing 
finite-size effects, I have chosen to simulate cubes with 
spins, Lx ranging from 10 to 80. For comparison, 
a few simulations have also been done for samples with 
anisotropic shapes. 

Critical phenomena are best studied numerically with 
non-local Monte Carlo methods. Recently, we proposed 
a cluster-flip algoijithm, based on the one-cluster flip al- 
gorithm of Wolff,c3 especially designed for jjiimulating 
spin systems with competing interactions.EZIo Below, 
this algorithm, which combines the cluster algorithm 
approaches, for systems with long-range ferromagnetic 
couplingscJ with that|Jor systems with nearest neigh- 
bor random couplings ,c3 is reformulated in order to take 
modified surface couplings into account. 

Let i be a lattice site characterized by the spin Si and 
already belonging to the cluster we iteratively build up. 
A nearest neighbor site j with spin Sj is added to the 
cluster with probability 



i (l-t- sign (s,sj)) (1 



cxp [-2J,/{kBT)]) 



(2) 



if sites i and j are both located in the surface layer, and 
with probability 



- (H-sign(siSj)) (1 



exp[-2J6/(fcsT)]) 



(3) 



otherwise. A lattice site k with spin s^. axial next-nearest 
neighbor to i is included with probability 



1 



(1 - sign (s.Sfe)) (1 - exp [-2 J,«;,/(fcBr)]) (4) 



During the simulation, layer-dependent quantities are 
computed. I discuss these quantities here only for case 
A, the corresponding quantities for case B are obtained 
by replacing z by y and by Ly. Of great interest 
when studying surface properties are the magnetization 
per layer 
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and the susceptibility per layer 
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The surface magnetization is then mi 
m(z = L^), whereas the response of the surface magne- 
tization to a surface field is xii = x{z = 1) = = ^z)- 
From the profiles m{z) and one also obtains the 

surface excess quantities 



E 

2 = 1 



(mf, - m{z)) 



and 



(8) 
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where rrifc and Xb sse the bulk magnetization and the 
bulk susceptibility. The response of the surface to a bulk 
field, xi, the energy, the specific heat, and the Binder 
cumulant have also been computed. 

Thermal averages are extracted by generating, after 
equilibration, 5 x 10^ clusters, and error bars are obtained 
by averaging over at least ten different realizations using 
different random numbers. 



III. SURFACES PERPENDICULAR TO THE 
AXIAL DIRECTION 



for case B if both spins are surface spins, and with prob- 
ability 

i (1 - sign {s,Sk)) (1 - exp [-2JbHb/ {kBT)]) (5) 

otherwise. The final cluster, which is flipped as a whole- 
presents two pecularities which are worth mentioning.c2l 
First, it contains spins of both signs, in variance with 
the Wolff cluster method, where all spins belonging to 
a cluster to be flipped have the same sign. Second, the 
cluster spins are not always connected by nearest neigh- 
bor bonds, which is again different from the traditional 
cluster flip method. This algorithm works well in the 
vicinity of the Lifshitz point, as demonstrated in the re- 
cent numerical study of the bulk Lifshitz point critical 
behaviour .EH 



Semi-inflnite ANNNI models with surfaces oriented 
perpendicular to the direction of competing interactions 
exhibit a crossover from Ising surface critical behaviour 
for Kf, = to Lifshitz point surface criticality for Kb = k^. 
For the purpose of studying this crossover, simulations 
were not only done for the value of Kb at the Lifshitz 
point, K^ = 0.27, but also at Kb = 0.15 and Kb = 0.24. 
The corresponding critical temperatures are listed in Ta- 
ble L Note that the same values of Kb haveJpeen used in 
a study of ANNNI bulk critical behaviouro 

Various values of the surface coupling strength, Jg, 
have been considered, with Jg ranging from to 3Jfc. The 
resulting surface phase diagrams are displayed in Figure 
2, together with the corresponding phase diagram of the 
semi-infinte Ising model, Kb = 0. The horizontal lines in- 
dicate the bulk critical temperatures which decrease with 
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increasing ki,. One also observes a shift of the location of 
the special transition point to lower values oir = Jg/Jb 
for larger values of the competition parameter Kb, in ac- 
cordance with the recent mean-field results.E3 This shift 
is mainly due to the decrease of the bulk critical temper- 
ature. The different surface transition lines in Figure 2 
merge for Jg S> Jb and form one line independent of the 
value of Kfc. This follows from the fact that for very strong 
surface couplings the surface effectively decouples from 
the bulk. At the Lifshitz point value Kb = 0.27, the phase 
diagram is also similar to the Ising case, with an ordinary 
transition, a surface transition, an extraordinary transi- 
tion, and a special transition point. This,is in agreement 
with the mean-field treatment Q£-|Ref.ca, but disagrees 
with the results obtained in Ref.L3. Based on our data, 
the special transition point of the semi-infinite ANNNI 
model at the Lifshitz point is located at r^^ — 1.15±0.05. 

Crossover phenomena showing up in the surface mag- 
netization may be analysed by plotting the effective ex- 
ponent 

Pi^eff^dXnmi/dhit (10) 

where t = (Tc — T)/Tc is the reduced temperature. On 
approaching Tc, t — > 0, /?i,e// becomes the critical ex- 
ponent Pi of the surface magnetization. In the following, 
data unaffected by finite-size effects are usually displayed, 
finite-size dependences-Jpeing circumvented by adjusting 
the size of the samplers 

Effective surface magnetization exponents obtained for 
Kb = 0, 0.15, and 0.24, with Jg — J\, are shown in Figure 
3. Clearly, for a fixed value of t, Pi^eff decreases with 
increasing Kt,. In the limit t — > 0, however, the effective 
exponents all tend asymptotically to the critical surface 
magnetization exponent f3i w 0.80 of ptiie semi-infinite 
Ising model at the ordinary transitionJ22 in accordance 
with the IsiBft character of the ANNNI bulk transition 
for Kb < k^.e3 The observed increase of the corrections 
to scaling with Kb reflects not only the closer proximity 
of the bulk Lifshitz point but also the reduced distance 
of r = Jg/Jb = 1 to the special transition point rsp{Kb)- 
This is illustrated in Figure 3 by the data obtained for 
Kb = 0.24 with Jg — 0.5Jb where corrections to scaling 
are greatly reduced compared to the case Jg = Jb 

Some of the results obtained for Kb = k^ , with r < r^p, 
are displayed in Figures 4 and 5. For this choice of the in- 
teractions, the surface undergoes an ordinary transition 
at the Lifshitz point critical temperature T^. Both the 
magnetization per layer, shown in Figure 4, and the sus- 
ceptibility per layer (not shown) vary non-monotonously 
close to T^. Whereas a non-monotonic behaviour of the 
layer susceptibility is also observed close to the ordinary 
transition in thick Ising films, the layer magnetization in- 
creases monotonicajly in Ising films from its surface value 
to the bulk value. t3 The non-monotonic behaviour of the 
layer magnetization close to the Lifshitz point ordinary 
transition may be explained in the following way. Inside 
the bulk the value of the magnetization mainly results 



from the balancing of competing influences in axial di- 
rection. Close to the surface, however, some negative 
contributions are missing, which then yields a tendency 
to stronger local ordering, resulting in a maximum of 
the layer magnetization near the surface. Similar non- 
monotonic profiles also show up in the semi-infinite crit=. 
ical Ising model in the presence of a weak surface field.a 
It should be noted that in our case the appearance of 
non-monotonic profiles depends strongly on the value of 
the competition parameter Kb. This behaviour is not ob- 
served for the studied values of Kb smaller than k^ . For 
i^b = the maximum in the magnetization profile shows 
up for all values of the surface coupling Jg leading to the 
ordinary transition. When the temperature is increased, 
the maxima both of the layer magnetization and of the 
layer susceptibilty are shifted towards the center of the 
system. 

The surface magnetizations and the corresponding ef- 
fective exponents obtained for two different values of the 
coupling ratio Jg/Jb are compared in Figure 5. As ex- 
pected, the amplitudes (see Figure 5a) and the correc- 
tions to scaling (see Figure 5b) differ. It is only close to 
Tc, for t < 0.07, that the effective exponents (and there- 
fore the corrections to scaling) become identical. Further- 
more, the effective exponents then vary almost linearly 
with temperature. A linear extrapolation yields the value 
Pi — 0.618 ±0.005 for the surface magnetization critical 
exponent at the Lifshitz point ordinary transition. This 
value is clearly smaller than the value Pi — 0.80 ± 0.01 
obtained ati-the ordinary transition of the semi- infinite 
Ising model.ES It is worth noting that mean-field approx- 
imation yields the common value Pi^^ = 1 for both cases. 

Estimates of Lifshitz point surface critical exponents at 
the ordinary transition are gathered in Table II. Hereby, 
the critical exponent 7^ = 0.84±0.05 of the susceptibilty 
Xi, i.e. of the response of the surface to a bulk field, can 
be obtained by analysing the corresponding effective ex- 
ponent in the same manner as discussed previously. The 
response of the surface to a surface field, xii, however, is 
finite at T^, with a cusp-like singularity when approach- 
ing the critical point from high temperatures. Instead 
of trying to extract the corresponding critical exponent 
'Yii directly from xii, it is easier to study the derivative 
xii = dxii/dT which in our case diverges on approach to 
the critical point with the power law x'u ~ (^^) 
This then yields the value -0.06 ± 0.02 for 7^^. 

Further Lifshitz point surface critical exponents are ob- 
tained by analysing the power law behaviour of the ex- 
cess quantities and Xs, see Eq. (8) and (9), close to 
the critical point: mg ~ and Xs ~ ■ Based on 
our data, the excess quantities critical exponents are es- 
timated in the present case to be P^ — —0.14 ± 0.04 and 
7^ = 1.69 ± 0.07. 

A closer inspection of Table II reveals that various 
scaling relations are fulfilled, thus demonstrating the 
reliability of our estimates. Indeed, inserting the Lif- 
shitz point bulk critical exponents 7^ — 1.36 ± 0.03 and 
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= 0.238 ± 0.0050 as well as i^f = 0.348,0 the follow- 
ing scaling relations are readily verified: 

Ps = 0t - (11) 

Is = it + y\ (12) 

27f-7n-7f (13) 

Lifshitz point surface quantities have also been stud- 
ied near the special transition point with J, = 1.15 J^, see 
above. At temperatures below T^^ magnetization profiles 
are extremely flat: only very close to the surface does one 
observe a slow increase of the layer magnetization which 
reaches a maximum in the layer immediately below the 
surface layer. The surface magnetization itself is slightly 
smaller than the bulk magnetization. Estimated values 
of Lifshitz point surface critical exponents at the special 
transition point are listed in Table II. As both the bulk 
and the surface become critical at the special transition 
point, the responses of the surface to a bulk field as well 
as to a surface field diverge. The corresponding critical 
exponents are estimated to be 7^ = 1.28 ± 0.08 and 
7ii sp — 0.76 ±0.05 where effective exponents have again 
been analysed. For the surface magnetization critical ex- 
ponent the value /3i,sp — 0.22 ± 0.02 is obtained. Error 
bars given for the special transition point critical expo- 
nents include the uncertainty in the localisation of this 
multicritical point. 

It is worth noting that the Lifshitz point surface criti- 
cal exponent at the special transition point given in ICfu 
ble II agree within the error bars with the cstimateal3'Q 
of the corresponding exponents for the semi-infinite Ising 
model. The change from Ising to Lifshitz point bulk crit- 
ical behaviour when k\, — s- seems to have only a mi- 
nor impact on the computed surface critical behaviour at 
the special transition point. However, this may not hold 
for the crossover exponent, which describes tlie merging 
of the surface transition line with this point ,u as differ- 
ent values are obtained in mean-field aperoximation for 
the Ising and for the Lifshitz point caseE3 In the present 
work, no attempt was made to determine the crossover 
exponent. 

IV. SURFACES PARALLEL TO THE AXIAL 
DIRECTION 

Two different approaches have been chosen for study- 
ing the critical behaviour of surfaces oriented parallel to 
the axial direction. In the first approach the value of the 
surface competition parameter Ks is set to zero so that 
only the strength of the surface nearest neighbor coupling 
Js is varied. In the second approach the bulk and surface 
competition parameters have the same value. Kg = Kf,, 
with Js ranging from to 6Jf,. Simulations done for the 



latter case may be viewed as a preparatory work for the 
study of the more complex multicritical points mentioned 
in Section II. In this work I mainly discuss Lifshitz point 
surface criticality, i.e. k;, = k^, at the ordinary transition 
and at the special transition point. 

First one has to note that the critical value of the cou- 
pling ratio = Jsl Jb, needed for the surface to get 
critical by itself, now depends both on the bulk and on 
the surface competition parameters. For example, at the 
bulk Lifshitz point r^^ — 1.75 ± 0.05 for = 0, whereas 
for Us — the special transition point is tentatively 

located at r£, « 4.3. This remarkable shift of the criti- 
sp 

cal coupling ratio to larger values in the latter case may 
be explained by a stronger weakening of the effectice sur- 
face couplings, due to the competing surface interactions, 
than of the effective bulk couplings. 

At the bulk Lifshitz point, the magnetization per layer 
increases monotonically from its surface value to its bulk 
value, as shown in Figure 6. This is in marked contrast 
to the behaviour encountered when the surface is per- 
pendicular to the axial direction, as discussed in the pre- 
vious Section, see Figure 4. For the surface orientation 
considered here, the layers are parallel to the direction of 
competing interactions, so that no axial contributions are 
missing close to the surface and no tendency to stronger 
local ordering is expected. 

Lifshitz point surface critical exponents at the ordinary 
transition, which have been computed in a similar way as 
for case A, are gathered in Table III. These values should 
be compared with the corresponding values of case A 
surfaces listed in Table II. Figure 7 shows effective expo- 
nents derived from the surface magnetizations obtained 
by simulating samples with different surface orientations. 
Clearly, both types of surfaces yield different asymptotic 
values Pi for the effective exponent, thus demonstrating 
that the values of surface critical exponents depend on 
the surface orientation. Different values are also obtained 
for the critical exponents j^^, (3^, and 7^". For 7^ the 
situation is not so clear as rather similar values follow 
from the analysis of effective exponents. One may note 
that mean-field theory yi£lds. for the same value for 
both surface orientationsHo 

Whereas the scaling relation (]l^) is fulfilled indepen- 
dently of the orientation of the surface, the scaling rela- 
tions ([1^) and ( p^ ) have to be modified for surfaces paral- 
lel to the axial direction. Indeed, the behaviour of excess 
quantities is governed close to a bulk critical point by the 
bulk correlation length along the direction perpendicular 
to the surface. The Lifshitz point being an anisotropic 
critical point characterized by two correlation lengths di- 
verging with different critical exponents, respectively 
has to be used for case B respectively A surfaces. This 
then yields for surfaces parallel to the axial direction the 
scaling relations 

Pt = Pt - (14) 

and 
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ls=lb+'^^- (15) 

Using = 0.709, a nice agreement with the numerical 
estimates for (3^ and is observed. 

I also studied the Lifshitz point surface critical be- 
haviour of surfaces oriented parallel to the axial direction 
close to the special transition point, choosing Ks = and 
Js = 1.75Jb, with Kb = . The estimates for the critical 
exponents Pi^p, Jisp: and sp have been included in 
Table III. Again, the quoted error bars take into account 
the sample averaging as well as the uncertainty in the 
location of the special transition point. Interestingly, the 
values of the critical exponents are very similar to the val- 
ues obtained for the other type of surfaces considered in 
Section III and quoted in Table II. This is a strong indi- 
cation that at the bulk Lifshitz point the surface critical 
behaviour at the special transition point may not depend 
on the surface orientation. 



V. CONCLUSION AND OUTLOOK 

Whereas surface criticality of models with only ferro- 
magnetic couplings has been the subject of numerous in- 
vestigations, studies of surface critical behaviour in sys- 
tems with competing interactions are very scarce. The 
aim of the present work is to fill some gaps by presenting 
the first Monte Carlo study of surface critical behaviour 
in semi-infinite three-dimensional ANNNI models, pay- 
ing special attention to the ordinary transition as well 
as to the special transitimipcoint, both at the bulk Lif- 
shitz point. Prior worksEJo exclusively focused on the 
mean-field treatment of this problem. 

Overall, the numerical results, which have been ob- 
tained by studying different surface orientations with re- 
spect to the axis of competing interactions, are in qualita- 
tive agreement w ith re sults obtained in the recent mean- 
field treatments.Ha The location of the special tran- 
sition point is a function of the bulk competition pa- 
rameter as well as of the surface competition parameter, 
when present. At the bulk Lifshitz point, the surface 
diagram is similar to that of the semi-infinite 3D Ising 
model, with a surface transition line ending in a special 
transition point-. This is in accordance with recent mean- 
field resultsjpo but disagrees with the early treatment 
of Gumbs.t3 Finally, at the bulk Lifshitz point, surface 
critical exponents at the ordinary transition depend on 
the orientation of the surface, as predicted by mean-field 
theory. 

The main results of our study are the Monte Carlo es- 
timates of surface critical exponents at the bulk Lifshitz 
point gathered in Tables II and III. This is the first time 
that values for these exponents have been determined 
numerically. Whereas at the ordinary transition surface 
critical exponents depend on the surface orientation, the 
computed values at the special transition point do not 
show a clear dependency on the orientation of the sur- 
face with respect to the axial direction. Furthermore, 



these values agree rather well with the values obtained 
in the Ising model at the special transition point. In 
order to clarify the situation at this multicritical point 
renormalization group calculations are desirable. How- 
ever, with regard to the immense technical difficulties 
already encountered in the field-theoretical coipoutation 
of bulk critical exponents at the Lifshitz point J25l3 this 
seems to be a formidable task not easily achievable in the 
near future. 

A host of open and interesting problems remains for 
future work. In the present study, surface criticality has 
only been investigated for values of the bulk competition 
parameter Kb < k^, thus neglecting the interesting case 
when the bulk undergoes a phase transition to a mod- 
ulated structure. Similarly, for surfaces parallel to the 
axial direction, the surface competition parameter Ks has 
been chosen in such a way that at the surface transition 
a uniformly ordered phase appears in the surface layer. 
However, for large values of Ks one will observe a surface 
transition to a floating incommensurate phase in the sur- 
face layer. Various interesting scenarios follow, one of the 
most intersting being a multicritical point where a sur- 
face transition line with a floating phase in the surface 
layer meets an ordinary transition line with a modulated 
bulk. 

A further interesting question concerns the scaling 
form at the Lifshitz point of spin-spin correlation func- 
tions along and perpendicular to the surface. A re- 
cent generalization of conformal invariajace to the strong 
anisotropic scaling at the Lifshitz pointc3 has yielded pre- 
dictions for the scaling form of bulk two-point correlators 
which were found to agree perfectly with Monte Carlo 
results.L3 Similarly, predictions for two-point correlation 
functions in semi-infinite systems at the Lifshitz point 
can be obtained in the framework of this theory and may 
then be compared to numerical data. Work along the 
sketched lines is planned for the future. 
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TABLE I. Bulk critical temperatures obtained for different 
values of the bulk competition parameter Kb- 





0.15 


0.24 


0.27 


ksTc/Jb 


4.131 ± 0.004 


3.865 ± 0.003 


3.7475 ± 0.0005 



TABLE II. Surface critical exponents at the bulk Lif- 
shitz point obtained for the 3D semi-infinite ANNNI model 
with the surface layer oriented perpendicular to the axial 
direction.!— 1 OT-MF: mean-field exponents at the ordinary 
transitionEJ, OT-MC: Monte Carlo values at the ordinary 
transition, SP-MC: Monte Carlo values at the special transi- 
tion point. The numbers in brackets give the estimated error 
in the last digit. 







7f 


7fi 


/3f 


7." 


OT-MF 


1 


1/2 


-1/4 


1/4 


5/4 


OT-MC 


0.618(5) 


0.84(5) 


-0.06(2) 


-0.14(4) 


1.69(7) 




Hl,sp 




lll.sp 






SP-MC 


0.22(2) 


1.28(8) 


0.76(5) 







TABLE HI. The same as in Table II, but now for the 3D 
semi-infinite ANNNI model with the surface layer oriented 
parallel to-fche axial direction. Mean-field values are extracted 
from RefEj. The mean-field values of thp-critical exponents 
f3^ and 7^^, which are not given in Ref.tj, follow from the 
assumption that the scaling relations ( |l^ ) and jl^ ) hold. 





/3f 


7f 


7fi 


I3s 


7." 


OT-MF 


1 


1/2 


-1/2 





3/2 


OT-MC 


0.687(5) 


0.82(4) 


-0.29(6) 


-0.46(3) 


1.98(8) 




Hl,sp 




7ii,sp 






SP-MC 


0.23(1) 


1.30(6) 


0.72(4) 
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FIG. 1. Cross sections of semi-infinite three-dimensional 
ANNNI models showing the two different types of surfaces 
studied in the present work: (a) surfaces perpendicular to the 
axis of competing interactions, (b) surfaces parallel to this 
a:xis. Jb and Js denote the nearest neighbor bulk and surface 
couplings, respectively, whereas the axial next-nearest neigh- 
bor interactions axe labeled by the bulk, Kb, and surface, Ks, 
competition parameters. Surface lattice sites are represented 
by filled points. 
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extraordinary 
transition 
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FIG. 2. Surface phase diagrams obtained for difi^erent val- 
ues of the bulk competition parameter Kb'. (circles), 0.15 
(squares), 0.24 (diamonds), and 0.27 (triangles). Computed 
critical temperatures are shown as symbols, the filled sym- 
bols indicate the locations of the different special transition 
points. The shift in the bulk critical temperatures (horizontal 
lines) as a function of Hb is obvious. Note that for the Lifshitz 
point value of Kb (ft^ = 0.27) a special transition point and a 
surface transition line are still encountered. 
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-a K^=0.15, J=J^ 
-^K,=0.24,J=J, 
-^K,=0.24,J=0.5J, 
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FIG. 3. Surface magnetization effective exponents as func- 
tion of the reduced temperature t for different interaction 
sets as indicated in the legend. Note that in the limit 
t — » all curves extrapolate to the common asymptotic value 
/3i « 0.80, thus demonstrating that the ordinary transitions 
encountered for these interaction sets all belong to the same 
Ising surface universality class. In order to avoid finite-size 
dependences Monte Carlo systems with up to 120 x 120 x 40 
spins have been simulated. Error bars result from averaging 
over at least ten different realizations. 
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0.69 



0.67 
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0.63 



FIG. 4. Layer magnetizations obtained in the Monte Carlo 
simulation of an ANNNI film containing 120 x 120 x 40 spins, 
with Kb = Kb and Js = 0.75Jb- The two data sets correspond 
to two difTerent temperatures: ksT/Jb = 3.55 (circles) and 
ksT/Jb ~ 3.60 (squares). Shown is the most interesting part 
where a non-monotonic behaviour as function of the layer 
index is observed. The error bars are far smaller than the 
sizes of the symbols. 




0.75 




0.65 



02. 




0.55 



FIG. 5. (a) Surface magnetization versus temperature and 
(b) effective exponent /?i,e// versus t for Kb = k^ . The data 
shown have been obtained for two values of the surface cou- 
plings: Js = (circles) and Ja = 0.75Jb (squares). A lin- 
ear extrapolation of the effective exponents close to t = 
yields the estimate /3i = 0.618 ± 0.005 for the surface mag- 
netization critical exponent at the Lifshitz point. In (a) the 
dashed line indicates the Lifshitz point critical temperature 
ksT^ / Jb = 3.7475. Only data not affected by finite-size ef- 
fects arc shown. Systems with up to 120 x 120 x 60 spins have 
been simulated. 




FIG. 6. Layer magnetizations computed for ANNNI sam- 
ples with surfaces parallel to the axial direction at differ- 
ent temperatures: ksT/Jb — 3.55 (circles), ksT/Jb = 3.60 
(squares) and ksT/Jb = 3.65 (triangles). Films with 
60 X 60 x 60 spins have been simulated, with Ks — Kb — Kb 
and Js = Jb- The layer magnetization increases monoton- 
ically with the distance to the surface. Error bars are far 
smaller than the sizes of the symbols. 



9 



0.75 



0.70 - 




^ 0.65 - 




0.60 - 



0.55 



0.0 



0.1 



0.2 



0.3 



t 



FIG. 7. Comparison of surface magnetization effective ex- 
ponents obtained at the bulk Lifshitz point from ANNNI sam- 
ples with different surface orientations: surface perpendicular 
to the axial direction (squares) or parallel to that axis (cir- 
cles). In the former case Js = 0.75J6, whereas in the latter 
case fvs f^b ~ 

and Js ~ Jb- Extrapolating these data 
to t = yields different values for the Lifshitz point surface 
magnetization critical exponent at the ordinary transition. In 
order to circumvent finite-size effects, different system sizes 
have been simulated. The largest system considered contained 
120 X 120 X 60 spins (squares) or 80 x 80 x 80 spins (circles), 
respectively. 
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